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We present preliminary results on a search for squarks and gluinos in proton-antiproton collisions with a center-of- 
mass energy of 1.96 TeV and based on about 2.0 fb~ 1 of data collected by the CDF detector in the Tevatron Run 
II. Events with multiple jets of hadrons and large missing transverse energy in the final state are studied within the 
framework of minimal supergravity (mSUGRA) and assuming R-parity conservation. The results are compared to 
Standard Model (SM) predictions, and limits on gluino and squark masses are extracted. A specific search for the 
supersymmetric partner of the bottom quark produced from gluino decays is carried out using a sample of events 
with missing transverse energy and two or more jets in the final state where at at least one of the two leading jets is 
b-tagged. Good agreement is found between data and SM predictions, and limits on gluino and sbottom masses are 
extracted. 

1. INCLUSIVE SEARCH FOR GLUINO AND SQUARK PRODUCTION 

Supersymmetry (SUSY) [T] is regarded as one of the most promising fundamental theories to describe physics at 
arbitrary high energy beyond the Standard Model. In mSUGRA [2] symmetry breaking is achieved via gravitational 
interactions and the vast SUSY parameter space is reduced to only five parameters: the common trilinear coupling 
at the GUT scale Aq; the sign of the Higgsino mixing parameter sign(fj,); the ratio of the Higgs vacuum expectation 
values tan/3; the common scalar and gaugino masses at the GUT scale mg and m 1 / 2 , respectively. At the Tevatron, 
the production of squarks (q) and gluinos (g), super-partners of quarks and gluons, constitutes one of the most 
promising channels because of the strong couplings involved. The cascade decay of gluinos and squarks into quarks 
and gluons will result in a final state consisting of several jets plus missing transverse energy coming from the 
neutralinos, which leave CDF undetected. 

The inclusive search [3] for q and g is performed using 2.0 /6" 1 of data collected by the CDF detector in Run 
II [4J. A mSUGRA scenario with Aq = 0, /i < and tan/3 = 5 is assumed. The gluino-squark mass plane is scanned 
via variation of the parameters m (0 - 500 GeV/c 2 ) and m 1/2 (50 - 200 GeV/c 2 ). The PYTHIA [5\ Monte Carlo 
program with the ISASUGRA implementation is used to predict the SUSY spectrum at the TeV scale. Light flavor q 
masses are considered degenerate, while 2-to-2 processes involving stop (t) and sbottom (6) production are excluded 
to avoid strong theoretical dependence on the mixing in the third generation. Signal samples are normalized to 
the next-to-leading-order (NLO) cross section calculated with PROSPINO [Bj. The SM backgrounds are dominated 
by QCD multi-jets processes where comes from partially reconstructed jets in the final state. Other sources 
of backgrounds are Z and W production in association with jets, top and diboson production where fx is mostly 
due to the presence of neutrinos or muons in the final state. The estimation of both signal and background yields 
relies on MC techniques. The background model is tested in different control regions defined to enhance a particular 
background contribution. A good agreement is always observed between data and SM prediction in all the control 
regions. 

For best sensitivity across the rrig — rrig plane, three different analyses are developed, requiring at least 2, 3 or 4 jets 
in the final state. The rejection of distinct SM backgrounds relies on specific cuts on kinematic variables including 
the transverse energy of the leading jets, the J^t of the event, and the angle between the Ifir and the direction of 
the jets. Different selection criteria are defined for each final state to enhance the signal over background separation 
in different regions of the mSUGRA phase space. A detailed study of the systematic uncertainties is carried out 
for each final state: the dominant systematic on the signal efficiencies and SM background yields comes from the 
3% uncertainty on the jet energy scale (JES). This translates into a 15% and 30% uncertainty on the signal and 
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background estimation, respectively. Other systematic uncertainties are related to the modeling of the initial- and 
final-state radiation (ISR/FSR) and the choice of Renormalization Scale and PDFs. Finally, a 6% uncertainty on the 
integrated luminosity is also included. The number of expected SM events is 48 ± 17 for the 4-jets analysis, 37 ± 12 
for 3-jets and 16 ± 5 for 2-jets. The observed events in 2 fb^ 1 of data are 45, 38 and 18 respectively. No significant 
deviation from the SM prediction is detected and the results are translated into exclusion limits on gluino and squark 
production. The limits are computed with a Baycsian approach at 95% confidence level (CL) and including all the 
systematics and their correlations in the calculation. Depending on gluino and squark masses, cross sections in the 
range 1 pb to 0.1 pb are excluded (see Fig. |T| left ) . In an mSUGRA scenario for which rrig ~ rriq this analysis excludes 
masses up to 392 GeV/c 2 . Squark masses below 280 GeV are excluded in any case (see Fig. [l] right). 




Figure 1: left: Exclusion limits on gluino and squark production cross section when m g = rriq. right: Observed (red region) and expected 
(dashed-line) 95% CL exclusion limits on gluino/squark production in the m g — m q plane. Regions excluded by previous experiments are 
shown, together with no-mSUGRA solution region. 



2. SEARCH FOR SBOTTOM PRODUCTION FROM GLUINO DECAY 

Under the hypothesis rri~ h < rrig < rriq, where rriq is the mass of the light flavor squarks, the gluino pair production 
pp — > gg is one of the dominant SUSY processes at Tevatron. In the framework of the Minimal Supersymmetric 
Extension of the Standard Model (MSSM), the assumption m t , m^o > m~ h > leads to a branching ratio of 100% 
for the b — > b\ + decay. In this configuration and with the conservation of the R-Parity, the sbottom production 
through gluino decays has a very distinctive signature with 4 6-jets and large missing transverse energy ($t) due to 
neutralinos in the final state. 

The present analysis [7] is based on 2.5/& -1 of data collected with the CDF detector during Run II. The event 
selection follows the strategy developed for the inclusive analysis requiring a minimum fJ^ per event and at least 
two jets in the final state. To enhance the heavy flavor content, at least one of the two leading jets is required to 
be 6-tagged. The spectrum of MSSM events is predicted with PYTHIA while PROSPINO is employed to calculate 
the NLO cross section. The gluino and the sbottom masses are varied from 240 GeV/c 2 to 400 GeV/c 2 , and from 
150 GeV/c 2 to 350 GeV/c 2 , respectively. The neutralino mass is fixed to = 60 GeV/c 2 and the squark mass 
to rriq — 500 GeV/c 2 . The estimation of the dominant background contributions from QCD multi-jet production 
and light flavor jets mis-tagged as 6-quark comes directly from data. The remaining SM backgrounds, including 
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the tt production, Z and W production in association with jets, single top and diboson production, are modeled 
using PYTHIA. The SM prediction is tested in background dominated control regions in which signal contribution is 
negligible. Good agreement between data and SM is always observed. Different sources of systematic uncertainties 
are studied. The primary contributions come from the JES (up to 25% on signal and background) and the mistags 
estimation (around 5%). 

Final kinematic cuts on the I^t of the event and the energies of the two leading jets are optimized to have the 
best sensitivity in two different areas of the SUSY phase space defined according to the value of Am = |m§ — mg|. 
In the small Am ~ 20 GeV/c 2 signal region 19 events are observed with a SM expectation of 22.0 ± 3.6. In the large 
Am ~ 80 GeV/c 2 signal region the expected SM yield is 22.7 ± 4.6 against 25 events observed in the data. Since 
no significant deviation from the SM prediction was observed, the results are translated into an exclusion limit on 
sbottom production process. Limit is computed with a Bayesian approach at 95% CL and including all the systematic 
uncertainties and their correlations in the calculation. Cross sections down to 40 fb are excluded. Results are also 
translated into the exclusion limit on gluino and squark masses shown in Figure [2j sbottom masses up to 320 GeV 
are excluded for mg\ u ~ 335 GeV. 




Figure 2: left: I$t distribution in the large Am signal region. Data are compared to expected SM and signal yields, right: Expected 
and Observed exclusion limits at 95% CL on — rrig plane. Limits from previous version of the analysis performed by CDF and D0 
during Run I and Run II[8j [9] are also shown. 
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